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thresholding  at  a fixed  percentile,  yields  primitives 
that  are  quite  effective  in  texture  discrimination. 


The  support  of  the  U.S.  Air  Force  Office  of  Scientific 
Research  under  Grant  AFOSR-77-3271  is  gratefully  acknowledged, 
as  is  the  help  of  Kathryn  Riley  in  preparing  this  paper. 

*Permanent  address:  Instituto  de  Pesquisas  Espaciais, 

C.  P.  515,  S.  J.  Campos-12200, 

S.  Paulo,  Brazil 


i ►«  i 'j  k-’  A .1  ■ i am  i it  I n tiu  *•  rr .. 


j?  REPORT  DOCUMENTATION  PAGE 

’ WSSHtET?  9-0  9 6 6/j1  “v 


2.  GOVT  ACCESSION  NO. 


4.  TljLZC«nd  SuI 


(£■  A COMPARISON  OF  SOME, ,§IMPLE METHODS  FOR  ) /S 

y.\  Extracting  texture  primitives  and  their  / / 

, ^EFFECTIVENESS  IN JEXTURE  ^DISCRIMINATION j 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  & PERIOD  COVERED 


Interim  ycc  * 


6.  performing  org.  report  number 


8.  CONTRACT  OR  GRANT  NUMBERf«J 


j Cj  ShvuanjWang,  Flavio  RtD 
Azriely Rosenfeld 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

University  of  Maryland 

Computer  Science  Center  1 ' 

College  Park,  Maryland  20742 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB , Washington,  DC  20332 

U.  MONITORING  AGENCY  NAME  & AODRESSf//  dlllerent  from  Controlling  Ollice) 


AFllSR— 77-3271  j 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 


61102F  2304/-A2 

l?.  a£RSBI.fiAT  E 

■ //  , Apr  a <879  j 

IV  NUMBER  OF  PAGES 

32 

IS.  SECURITY  CLASS,  (ol  I hi,  report) 


UNCLASSIFIED 

15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


*7.  DISTRIBUTION  ST.  IENT  (of  u * abstract  entered  in  Block  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  • fES 


I 19.  KEY  WOROS  (Continue  on  reverse  side  it  necessary  and  Identify  by  block  number) 


1 20  ABST^VC  T ( Continue  on  reverse  side  It  necessary  and  Identify  by  block  number) 


Three  simple  methods  of  extracting  texture  primitives  are  compared.  It 
appears  that  the  simplest  of  these,  thresholding  at  a fixed  percentile, 
yields  primitives  that  are  quite  effective  in  texture  discrimination.^ 

N 


DD  ,^nrm7,  1473 


^ 2 C / / 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  Of  TmiS  PAGE  . Hli^i  Mum 


1 . Introduction 

One  can  view  textures  as  composed  of  "primitives"  (connected 
regions  satisfying  certain  properties)  placed  in  a certain 
spatial  arrangement.  In  this  view,  to  describe  a texture  one 
needs  to  describe  both  the  primitives  and  the  placement  rules. 

For  the  primitives  we  can  construct  a property  list,  while  for 
the  placement  rules  we  can  select  a set  of  spatial  relations  and 
list  which  pairs  (triples,  quadruples,...)  of  primitives  satisfy 
these  relations.  This  model  has  been  used  by  many  researchers 
initially  to  synthesize  textures  [1]  , [2]  and,  more  recently,  in 
texture  discrimination  [3],  [4],  [5]. 

In  Tsuju  and  Tomita  [3],  the  problem  was  to  segment  a picture 
consisting  of  regions  of  different  textures.  The  primitives 
("atomic  regions")  of  which  the  textures  were  composed  were  de- 
fined to  be  connected  sets  of  points  with  almost  the  same  gray 
level.  Once  the  atomic  regions  were  extracted,  for  each  one  a 
set  of  properties  was  measured  (shape,  size,  position,  color, 
and  average  gray  level).  Histograms  of  the  property  values  were 
constructed  and  the  textures  were  identified  as  peaks  in  these 
histograms.  Having  identified  the  textures,  the  picture  was 
segmented  by  labeling  the  primitives  with  the  names  of  the  tex- 
tures they  belonged  to. 

The  pictures  to  which  Tsuji  and  Tomita  applied  their  methods 


consisted  of  artificially  constructed  objects  covered  by  almost 
regular  subpatterns.  An  extension  of  the  method  to  more  realistic 


scenes  was  proposed  by  Zucker  et  al.  [4].  For  "real-world" 
pictures  it  was  argued  that  the  notion  of  "atomic  region"  was 
not  well  defined,  or  better,  that  its  definition  did  not  lead 
easily  to  its  extraction.  Therefore,  instead  of  segmenting  the 
picture  to  extract  the  primitives,  local  property  values  were 
directly  obtained  through  a set  of  "spot  detectors"  having  a 
range  of  sizes.  What  was  measured  was  the  output  intensity  of 
the  spot  detector  of  a given  size  at  all  points  in  the  picture. 
To  prevent  "spurious"  spot  detections,  a spot  value  was  ignored 
if  a larger  spot  value  occurred  in  its  receptive  field.  In  this 
way  one  can  construct  a histogram  of  the  outputs  of  the  spot 
detector.  As  in  [3],  the  textures  were  identified  as  peaks  in 
this  histogram.  Assuming  that  we  have  two  textures  in  the 
picture,  they  can  be  retrieved  by  thresholding  the  histogram, 
e.g.  at  the  lowest  point  between  the  two  peaks.  The  spot  size 
was  selected  as  the  one  that  gave  the  most  strongly  bimodal 
histogram.  Although  the  method  based  on  spot  detection  gives 
less  detailed  information  about  the  primitives  than  the  one 
proposed  in  [3] , it  is  sufficient  for  discriminating  some  real- 
world  textures. 

In  both  papers  ([3]  and  [4])  no  attempt  is  made  to  use  the 


spatial  arrangement  of  the  primitives  to  discriminate  the  textures 
(it  was  not  necessary).  In  the  work  of  Maleson  et  al.  [5], 
however,  both  sources  of  information  are  suggested  for  texture 
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discrimination.  To  simplify  the  description  of  the  primitives, 
they  are  restricted  to  be  clusters  of  pixels  having  "simple" 
shapes.  As  primitive  properties,  it  is  suggested  that  average 
intensity  (gray  level),  eccentricity,  axial  orientation,  and 
size  constitute  a sufficient  set  for  most  textures.  The  place- 
ment of  primitives  is  described  by  a restricted  set  of  spatial 
relations  (such  as  collinearity ) between  primitives  in  the  same 
class  (i.e.,  having  the  same  property  values  for  a subset  of 
the  properties) . 

In  our  work  we  propose  three  related  methods  for  extracting 
texture  primitives.  These  methods  are  compared  in  connection  with 
the  discrimination  of  a few  texture  samples,  including  four 
Brodatz  [6]  textures  and  three  LANDSAT  geological  terrain  tex- 
tures. In  the  methods  proposed  the  primitives  are  not  forced  to 
conform  to  pre-specif ied  shapes  (as  they  were  in  (5]).  This  makes 
possible  the  use  of  shape  attributes  not  present  in  [5]  that  may 
be  useful  in  texture  discrimination.  Also  we  have  not  used 
spatial  relationships  between  the  primitives  to  characterize  the 
textures . 
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. Three  ways  of  extracting  texture  primitives 

The  basic  idea  of  all  the  methods  is  to  segment  the  texture 
into  primitives  and  background  by  thresholding  of  the  gray  level. 

We  will  define  the  primitives  as  being  regions  formed  by  con- 
nected sets  of  pixels  whose  gray  level  values  are  above  (or 
below)  the  threshold. 

Of  course,  by  thresholding  we  should  not  expect,  in  general, 
to  obtain  exactly  the  perceptual  texture  primitives.  But  it  is 
hoped  that  the  regions  obtained  are  sufficiently  related  to  the 
"actual"  texture  primitives  that  the  properties  measured  for  these 
regions  reflect  property  values  of  the  "actual"  primitives. 

We  can  predict  that  this  method  will  be  computationally  effi- 
cient, since  connected  components  of  the  binary  image  produced 
by  thresholding  can  be  identified  in  a single  pass  through  the 
picture.  Also  in  the  same  pass  much  information  about  the  primi- 
tives can  be  gathered  (area,  perimeter,  moments...). 

The  main  problem  is  how  to  select  an  adequate  threshold.  If 
the  histogram  of  the  gray  level  values  is  bi-modal,  a natural 
candidate  for  a good  threshold  is  the  valley  between  the  two 
peaks.  To  make  easier  the  computation  of  the  primitive  properties 
we  would  like  the  primitives  to  be  isolated  regions.  If  the 
threshold  that  separates  the  two  peaks  is  near  the  50th  percentile, 
most  likely  we  will  not  get  isolated  regions.  Also  if  the  histo- 
gram is  not  bi-modal,  the  approach  is  not  applicable.  Therefore, 


instead  of  making  the  threshold  depend  on  the  histogram,  in 
the  first  method  we  chose  the  threshold  such  that  we  have  a 


fixed  percentage  of  the  pixels  with  above  threshold  gray  level. 
The  percentage  must  be  low  enough  so  we  can  expect  to  get  a 
set  of  isolated  primitives,  but  sufficiently  high  to  make  the 
primitive  properties  meaningful.  We  have  chosen  25%,  but  other 
values  near  25%  would  be  good  as  well. 

For  each  texture  we  have  two  possible  thresholds  depending 
on  whether  we  choose  to  consider  dark  or  bright  primitives. 
Actually  we  may  use  both  types  of  primitives  for  texture  discri- 
mination. The  thresholding  of  the  texture  can  be  preceded  by 
some  pre-processing  operations,  such  as  noise  cleaning  or  even 
blurring.  Gray  level  transformations,  however,  such  as  histo- 
gram flattening,  do  not  significantly  alter  the  final  output. 

The  second  method  of  extracting  texture  primitives  is  based 
on  the  histogram  peak  sharpening  algorithm  of  Peleg  [7],  This 
algorithm,  which  can  be  thought  of  as  a one-dimensional  cluster- 
ing algorithm,  iteratively  alters  the  histogram  until  it  consists 
of  a few  spikes.  The  picture  is  segmented  into  the  few  gray 
levels  corresponding  to  the  spikes.  The  number  of  spikes 
(clusters)  of  the  modified  histogram  is  a function  of  the  number 
of  iterations  and  the  "neighborhood"  size  used  in  the  clustering 
process.  These  properties  were  used  to  produce  a set  of  spikes 
that  sums  up  as  nearly  as  possible  to  25%  of  the  picture 
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points.  In  this  way  we  take  into  account  characteristics  of 
the  histogram,  which  was  not  done  in  simple  thresholding. 

The  third  and  last  method  (" Superslice")  was  developed  by 
Milgram  [8].  In  it,  the  image  is  not  segmented  by  using  a 
single  threshold;  rather,  a collection  of  thresholds  is  used 
to  extract  regions  which  are  compared  to  a thinned  edge  map 
of  the  picture.  The  regions  are  accepted  or  rejected  based 
on  the  coincidence  of  the  boundary  points  of  the  regions  with 
edge  points.  Surviving  competing  regions  are  compared  and  only 
the  best  match  (with  the  edge  map)  is  retained.  This  method 
not  only  uses  histogram  information  but  also  complements  it 

; | 

with  local  information  derived  from  the  picture.  In  our  case 
we  used  thresholds  at  the  percentiles  15%,  25%,  and  35%. 

For  each  method  the  output  is  a binary  picture.  The  fixed 
threshold  method  is  the  computationally  cheapest,  while  Super- 
slice is  the  most  demanding.  The  last  two  algorithms  (iterative 
histogram  modification  and  Superslice)  are  described  briefly  in 
Appendix  A. 

In  order  to  get  meaningful  primitive  statistics,  we  use 
128x128  pictures.  The  textures  shown  in  Figures  la,  2a,  3a,  and 
4a  were  taken  from  Brodat  [6]  and  show  respectively  wool,  raffia, 
sand,  and  grass.  The  results  of  applying  thresholding  to  the 
noise  cleaned  pictures  are  shown  in  Figures  lc  to  12c.  Figures 
16a,  17a,  and  18a  show  three  different  terrain  types;  Mississipian 
limestone  and  shale,  Pennsylvanian  sandstone  and  shale,  and 


Lower  Pennsylvanian  shale.  These  pictures  were  noise  cleaned 
and  then  subjected  to  thresholding  (bright  25%  and  dark  25%) 
and  to  Superslice.  The  results  of  these  operations  are  shown 
in  Pictures  13b,  14b,  15b  (bright  25%) , 16b,  17b,  18b  (dark 
25%)  , 19c,  20c,  and  21c  (Superslice) . 


3 . Primitive  properties 


Once  we  have  extracted  the  texture  primitives,  the  next 
step  is  to  measure  properties  that  may  characterize  them.  We 
have  chosen  a set  of  shape  attributes  (area,  perimeter,  compact- 
ness, eccentricity,  direction)  together  with  average  gray  level. 

The  area  is  defined  simply  as  the  number  of  pixels  in  the 
primitive,  while  the  perimeter  is  its  number  of  boundary  pixels. 
Compactness  is  the  ratio  of  the  square  of  the  perimeter  to  the 
area.  Eccentricity  is  defined  as  the  ratio  of  the  major  to  the 
minor  axis  of  inertia  (the  detailed  formulas  are  given  in 
Appendix  B) . Finally,  direction  is  the  angle  that  the  major 
axis  makes  with  a fixed  orientation  (vertical  in  the  picture) . 

The  property  values  of  a given  region  serve  to  describe 
that  particular  region.  The  distribution  of  the  values  for  all 
regions  serves  to  characterize  the  texture;  it  can  be  thought  of 
as  a statistical  description  of  the  texture.  The  property  value 
distribution  is  given  by  the  property  value  histogram.  In  com- 
puting this  histogram  we  do  not  take  into  account  border  regions, 
i.e.,  regions  that  are  adjacent  to  a border  of  the  picture, 
becuase  such  regions  would  falsify  the  property  value  distribu- 
tion. * 

To  discriminate  the  textures,  we  use  features  derived  from 
the  histogram.  We  could,  for  instance,  take  as  features  the 
values  of  the  histogram  for  each  interval.  This  would  completely 
describe  the  histogram  but  would  give  too  large  a number  of 


samples  chosen  from  two  sets  of  textures.  The  first  consisted 
of  four  Brodatz  textures  [6]:  wool,  raffia,  sand,  and  grass; 
the  second  consisted  of  three  different  geological  terrains  from 
a LANDSAT  image.  In  each  set  there  were  four  samples  of  each 
texture . 

A feature  is  said  to  "separate"  two  given  textures  if  the 
ranges  of  feature  values  for  the  two  textures  are  non-overlapping 
intervals.  The  "separability"  of  two  textures  was  measured  in 
terms  of  the  number  of  features  that  separated  the  two  textures. 

For  the  Brodatz  textures  all  three  methods  described  in 
Section  2 were  used,  while  for  the  terrain  textures  only  (dark 
and  bright)  25%  thresholding  and  Superslice  were  used.  Figures 
21,  22,  23  display  the  feature  values  obtained  for  the  Brodatz 
texture  samples  (A,  B,  C,  D are  wool;  E,  F,  G,  H are  raffia; 

I,  J,  K,  L are  sand;  M,  N,  0,  P are  grass)  . (The  textures  de- 
picted in  Figures  la,  2a,  3a,  4a  are  A,  E,  I,  and  M,  respectively.) 
Tables  1,  2,  and  3 summarize  the  effectiveness  of  the  features 
for  separating  the  textures.  From  these  tables  we  see  that  no 
feauture  alone  separates  all  texture  pairs.  However,  two  features 
used  independently  are  sufficient  to  separate  all  pairs  of  tex- 
tures. 

A comparison  of  the  methods  can  be  made  if  we  look  at  the 
features  that  are  common  to  every  method.  This  information  is 


features.  We  decided,  instead,  to  use  the  mean  and  standard 
deviation  as  descriptors  for  the  distribution  of  the  property 
values.  In  the  computation  of  the  mean  and  standard  deviation 
very  small  regions  and  very  large  ones  were  not  considered. 

For  very  small  regions  most  of  the  properties  are  not  interest- 
ing; therefore  every  region  having  less  than  10  pixels  was  ig- 
nored. For  large  regions  we  used  a more  complex  rule  (it  is 
described  in  Appendix  C) . They  were  ruled  out  on  the  basis  that 
they  were  "atypical"  regions. 

The  effect  of  "erasing"  these  regions  from  the  binary  picture 
is  shown  in  Figures  Id,  2d,  3d,  4d,  5d,  6d,  7d,  8d,  9d,  lOd,  lid, 
12d,  13c,  14c,  15c,  16c,  17c,  18c,  19d,  20d,  21d.  (Actually 
the  regions  were  not  erased  but  just  ignored.) 

The  only  feature  derived  from  direction  was  standard  devia- 
tion. The  mean  was  not  used  because  it  would  vary  if  we  rotated 
the  texture.  Also  used  as  features  were  the  number  of  small 
regions  deleted  and  the  number  of  regions  used  in  the  computation 
of  the  mean  and  standard  deviation.  Thus  we  obtained  13  features; 
certainly  not  all  of  them  are  independent  of  the  others. 

For  the  Superslice  method,  only  ten  features  can  be  obtained. 
Since  Figures  18a,  19a,  20a,  21a  are  pseudo-images,  we  lose 
three  features:  the  number  of  small  regions  deleted,  the  mean 
gray  level,  and  the  standard  deviation  of  gray  level. 
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contained  in  the  upper  half  of  Table  7.  From  there  it  seems 
they  are  roughly  equivalent  with  a slight  advantage  to  simple 
thresholding. 

For  the  geological  terrain  textures  only  thresholding  and 
Superslice  were  used.  In  the  case  of  thresholding,  both  the 
bright  25%  and  the  dark  25%  were  used.  The  feature  values  are 
displayed  in  Figures  24,  25,  and  26  (A,B,C,D  are  Mississipian ; 

E,  F,  G,  H are  Pennsylvanian;  I,  J,  K,  L are  Lower  Pennsylvanian) . 
The  textures  depicted  in  Figures  12a,  13a,  and  14a  are  samples 
A,  E,  and  I.  The  effectiveness  of  the  features  is  summarized 
in  Tables  4,  5,  and  6 (for  Figures  24,  25,  and  26). 

In  this  set  of  textures  too  the  results  are  better  for 
simple  thresholding.  The  results  differ  from  those  obtained 
for  the  Brodatz  textures,  in  that  the  best  features  were  not 
the  same  ones.  This  should  case  no  surprise,  as  the  two  sets 
of  textures  are  quite  different.  Also  we  see  that  just  one 
feature  (mean  gray  level)  is  sufficient  to  separate  all  textures 


(Table  4) . 


5 . Discussion 

From  the  results  obtained,  it  seems  that  the  methods  sug- 
gested for  primitive  extraction  may  be  succesfully  employed 
for  texture  classification.  It  was  surprising  (but  pleasing) 
that  the  simplest  method  performed  better  for  both  sets  of 
textures . 

Although  all  texture  pairs  could  be  separated  with  two 
features  used  independently  (Brodatz)  or  one  feature  (terrain) , 
this  does  not  mean  that  for  larger  sample  sizes  we  would  have 
the  same  results.  However,  it  gives  us  an  indication  of  the 
usefulness  of  the  method.  Note  that  we  did  not  need  to  use 
pairs  or  triples  of  features  concurrently. 

Thresholding  and  iterative  histogram  modification  are  similar 
in  that  they  do  not  use  local  information  from  the  texture. 
Superslice  does  use  such  information  in  the  form  of  an  edge  map. 
The  binary  pictures  obtained  using  Supeislice  look  (subjectively) 
more  closely  related  to  the  original  textures  than  those  obtained 
using  the  thresholding  and  histogram  modification  techniques. 

In  spite  of  this,  Superslice  did  not  perform  as  well  as  regards 
texture  classification.  This  may  be  because  the  sample  size  was 
too  small,  or  because  of  the  dependence  of  the  features  (a  method 
performing  well  with  one  feature  will  also  perform  well  with  a 
correlated  one) . Further  studies  are  necessary  to  clarify  this 


matter . 
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Appendix  A 

Algorithms  for  Iterative  Histogram  Modification  and  Superslice 
Algorithm  A.l 

Iterative  Histogram  Modification  [7] 

Input:  Histogram. 

Let  be  the  number  of  pixels  having  gray  level  i.  For  each 
histogram  in  i,  the  neighboring  2r  (an  input  parameter)  bins 

i±j  on  each  side  of  i (j=l,2,...r)  are  examined.  If  is  greater 
than  the  average  A of  B^+1» . . . ,B^+r  (and  similarly  on  the  other 
side  of  i),  we  compute  the  ratio  X = (Bi~A)/B^,  which  specifies 
the  fraction  of  pixels  whose  gray  levels  will  be  shifted  towards 
i.  Then  the  following  gray-level  changes  are  executed: 

B^+r  • X from  i+r  to  i+r-1; 

B^+r_1  • X from  i+r-1  to  i+r-2;... 

B. • X from  i+1  to  i. 
l+l 

The  entire  process  is  iterated. 

Algorithm  A. 2 
Superslice  [8] : 

Input:  Gray  level  range  for  thresholding,  thinned  edge  picture. 

1.  For  each  gray  level  in  the  range: 
a)  threshold  the  image; 

label  all  connected  regions  of  above -threshold  points; 


b) 


c)  for  each  connected  region: 


i.  compute  the  percentage  of  border  points  which 
coincide  with  significant  thinned  edge  points; 

ii.  compute  the  contrast  of  the  region  with  the  back- 
ground ; 

iii.  classify  the  region  as  object/non-object  based  on 
the  size,  edge  match  and  contrast. 

2.  Construct  the  canonical  tree  for  the  set  of  object  regions 
based  on  containment. 

3.  Prune  the  containment  tree  by  eliminating  adjacent  nodes 
which  are  too  similar. 


Eccentricity  and  direction  can  be  computed  from  the  second 

order  moments  of  inertia  I , I and  I . Given  a set  of 

xx  yy  xy 


points  { * * ' * ^ f°r®ing  a region,  the  moments  are 
defined  to  be: 


Xxx  = 

,2 

ryy  = 

^x!2 

rxy 

.'ImiXjy^ 

where  x| 

and  y| 

are  the  coordinates  of 

point 

(xi'V 

relative  to  the 

center 

of  mass  of  the  region. 

The  angle  <p  which  the  larger  axis  makes  with  a horizontal 
line  is  given  by 


The  eccentricity  is  given  by 


Appendix  C 

Rule  for  Deleting  Large  Primitives 


*' 

The  area  range  of  the  region  for  a given  texture  sample  is 
divided  into  10  intervals;  each  interval  corresponds  to  an 
entry  in  the  histogram.  (In  this  histogram,  the  very  small 
primitives  and  the  primitives  touching  the  border  are  not  in- 
cluded.) The  upper  half  of  the  histogram  (from  the  sixth  to 
the  tenth  interval)  is  investigated  for  possible  deletions.  If 
we  have  more  than  three  regions  in  this  half,  no  deletion  is 
performed.  If  we  have  three,  every  region  (in  this  half) 
larger  than  400  pixels  is  deleted;  in  case  of  two  or  one,  the 
threshold  is  300  pixels. 

These  thresholds  (400,300)  should  not  be  thought  of  as 
absolute  numbers  but  as  percentages  of  the  total  window  area 
of  the  sample.  Considering  that  each  window  has  128x128  pixels, 
400  and  300  are  approximately  2.5%  and  1.9%  of  the  total  window 
size . 
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Table  1.  Effectiveness  of  the  features  for  the  Brodatz 
textures,  when  the  primitives  are  extracted 
by  thresholding  (25%).  A "Y"  entry  means 
that  the  ranges  of  the  feature  values  (cor- 
responding to  the  two  textures)  do  not  overlap 


I 


Table  3.  Analogous  to  Table  1 for  Superslice. 
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Table  2.  Analogous  to  Table  1 for  the  spike  program. 
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WITH  "Y" 

5 

2 

3 

10 

Table  4.  Effectiveness  of  the  features  in  dis- 
criminating pairs  of  terrain  textures. 
A "Y"  entry  means  that  we  have  non- 
overlapping ranges  of  feature  values. 
The  primitives  were  extracted  by 
thresholding  (bright  25%) . 
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Table  6.  Analogous  to  Table  4 but 


extracting  the  primitives 
using  Superslice. 


l.(a)  Wool;  (b)  complemented  wool; 
(c)  thresholded  (25%);  (d)  after 

the  deletion  of  small,  large, 
and  boundary  elements 
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Analogous  to  Figure  5 
for  sand 
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(a)  Wool  after  applying  Superslice 
the  lighter  regions  represent  the 
best  matches;  (b)  same  as  (a)  but 
complemented;  (c)  after  retaining 
the  best  matches;  (d)  after  delet- 
ing small,  large,  and  boundary 
regions 


13.  (a)  Complemented  Mississipian 
limestone  and  shale 

(b)  Thresholded  (25%  BRIGHT) 

(c)  After  deletion  of  small, 
large,  and  boundary  regions 


14.  Analogous  to  Figure  13,  for 

Pennsylvanian  sandstone  and  shale 


15.  Analogous  to  Figure  13,  for 
Lower  Pennsylvanian  shale 


16.  Analogous  to  Figure  13 
25%  DARK 


7 . Analogous  to  Figure  14 
25%  DARK 


, using 


8.  Analogous  to  Figure  15 
25%  DARK 


b 
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21.  Analogous  to  Figure  19,  for 
Lower  Pennsylvanian  shale 


